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We present an analytical model for the time-developing turbulent boundary layer (TD-TBL) over 
a flat plate. The model provides explicit formulae for the temporal behavior of the wall-shear stress 
and both the temporal and spatial distributions of the mean streamwise velocity, the turbulence 
kinetic energy and Reynolds shear stress. The resulting profiles are in good agreement with the DNS 
results of spatially-developing turbulent boundary layers at momentum thickness Reynolds number 
equal to 1430 and 2900 [1-3]. Our analytical model is, to the best of our knowledge, the first of its 
kind for TD-TBL. 



Formulation of the problem. In this paper we derive 
an analytical model for the time-developing turbulent 
boundary layer (TD-TBL) over a flat plate. We consider 
a flat plate (z = 0) submerged in an incompressible 
viscous fluid at rest for time t < 0. At time t = the 
fluid moves as a whole in the x-direction with velocity 
Voq. This motion creates a boundary layer near the 
plate. We assume that this boundary layer is turbu- 
lent. Note that in the case of TD-TBL all statistical 
characteristics of the flow depend only on the time t 
and distance z from the wall [e.g. the mean streamwise 
velocity is V(z,t). In contrast, in a spatially-developing 
TBL (SD-TBL) the statistical characteristics depend on 
two spacial directions, z and x, and the mean velocity 
has also V z component normal to the wall. Nevertheless 
in the limit of large Reynolds number both the TD-TBL 
and SD-TBL become asymptotically equivalent (with 
the replacement x Voot) [4]. Therefore it is reasonable 
to first consider the simpler simple case of TD-TBL. 

Definitions and Model Equations. We start from the 
Navier-Stokes equations for an incompressible fluid. The 
velocity U(r,t) is decomposed into the sum of its mean 
value V(z,t) = (U(r,t)) and a fluctuating part u(r,t), 
U(r, t) = xV(z, t)+u(r, t) Here x is the unit vector in in- 
direction, r = {x , y , z} is three dimensional coordinate, 
and (. . .) denotes averaging in time and in the span-wise 
direction y. 

The three main quantities in the model are the mean 
shear S(z,t), the tangential Reynolds stress r(z,t) and 
the turbulence kinetic energy per unit mass K(z,t), de- 
fined as: 



S(z,t) = ^ , T (z,t) = -(u x u z ) , K(z,t) = i (\u\ 2 ) . 

(1) 

The mean momentum equation ( e.g. Eq. (4.12) in [5]) 



after integration over z has the form: 



vS(z,t)+r(z,t) = r*(i) + 



d_ 

di 



V(z',t)dz' . (2) 



Here v is the kinematic viscosity and the right hand side 
(RHS) is momentum flux toward the wall. The integra- 
tion constant r*(t) = vS(0,t) is the wall shear stress. 

The turbulence kinetic energy conservation equation 
for a TD-TBL (see e.g. Eq. (5.132) in [5]) can be written 
as: 

dK(z,t)/dt + S(z,t) + V -T(z,t) = r{z,t) S(z,t) . (3) 

The RHS of this equation represents the production of 
turbulence kinetic energy by the mean shear. The two 
terms in the left hand side (LHS) : the rate of energy dis- 
sipation, £ , and spatial energy flux, T, require modeling 
via S, t and K . 

Equations (2) and (3) for S and K are exact. In or- 
der to solve them, we need to add a third equation for 
t(z, t) and model both £(z, t) and T(z, t). It is reasonable 
to assume that in the log-layer (region where the von- 
Karman log-law holds) the shear and normal Reynolds 
stresses have the same z and t dependence, thus their 
ratio is constant: r(z,t) / K(z,t) = c? . 

As suggested by [5] we write the rate of turbulence 
kinetic energy dissipation as £{z,t) = b[K(z,t)] 3 ^ 2 j z, 
where b is a positive constant and the spatial energy 
flux in the z direction as T(z,t) = —D(z,t) dK(z,t)/dz, 
where D(z,t) = dz yj K(z, t) is the turbulence diffusivity 
and d is a positive coefficient. 

Finally, we summarize the equations of the present 
model [Minimalist Model (Min-Model)] for TD-TBL as 



vS(z,t)+r(z,t) 
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d b^K{z~tj 
dt z 

T(z,t) 



r*(t) + J dz'—V(z',t), (4a) 

d d 
d—z^jK{z,t)— K(z,t) 
oz ozi 



t(z, t)S(z, t) 
c 2 K(z,t) . 



(4b) 
(4c) 
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The boundary conditions for Eqs. (4) at the wall (z = 0) 
and at the edge of the TBL [z = Z(t)) are: 

at the wall z = : V{Q,t) = 0, K(0,t)=0, (5a) 
at z = Z(t): V(Z, t) = Voo, K(Z(t),t) =0 . (5b) 

The numerical values of the model coefficients b and c 
are prescribed to be 0.34 and 0.53 respectively to obtain 
results in agreement with experiments ( e.g. [5]), and 
with DNS data (see Fig. 3 in [6]). The third coefficient 
d is fixed equal to 0.07 in order to match the role of the 
turbulent diffusivity in fully-developed turbulent channel 
flows [7]. 

We normalize the variables of the TD-TBL in "wall 
units" using the friction velocity u*(t), viscous length- 
scale and viscous time-scale t*(t) defined as 



(6) 



and denote the normalized variables V + = V/u«, K + = 
K/ul, z + = z/£*, t + = tit*. The friction velocity 
Reynolds number is defined as Re T (t) = u*Z(t)/v = 
Z(t)/e* = Z + {t), which is the width of TD-TBL in wall 
units. It should be noted that a similar normalization is 
usually employed for the spatially developing TBL (SD- 
TBL) except that the dependent variables in that case are 
functions of the streamwisc x location instead of time. 

The main advantage of the above normalization is 
that in wall units the stationary SD-TBL demonstrate 
universal behavior. We assume that the time- developing 
TBL exhibits similar universality in the limit of large 
times, as will be clarified later. We will also show that 
the time-developing TBL has properties very similar to 
that of the stationary TBL and its subregions can be 
classified in the same manner. 

"Logarithmic accuracy" for "asymptotically large 
times" For "asymptotically large times" we mean that 
lnt + ^> 1 . We assume that in this limit there is a re- 
gion in the TBL called the log-layer (with z + larger than 
upper boundary of the buffer sub-layer zi f ), where the 
mean streamwise velocity profile is described by the von- 
Karman law: 



V + {z + ) = K- L \nz+ + B, z + > z, 



buf 



50 



(7) 



where n w 0.41 is the von-Karman constant, and B w 
5.2. We define the edge of the TD-TBL, Z(t), as the 
normal distance from the wall where V = V^». Thus, the 
von-Karman law at Z + {t) becomes: V£, = V + (Z + ) = 
K~ 1 hiZ + + B. Moreover, in the limit \nt + 3> 1 the 
width of TBL Z{t) is large enough such that lnZ + 3> 
kB ~ 1, thus, lnRe r = lnZ + 3> 1. In the next section 
where we solve Eq. (4), this inequality will allow us to 
neglect terms of order of unity with respect to terms of 
order In Z + . We thus denote the accuracy of our results 
as the "logarithmic accuracy". For example, with the 
logarithmic accuracy, Eq. (7) for becomes nV^ — 
lnZ+. 



Self-consistent factorized solution. In order to re- 
duce the system of partial differential Eqs. (4) for the 
unknown functions S(z,t) and K(z,t) of two variables 
z and t to a system of ordinary differential equations it 
is convenient to introduce the following non-dimensional 
functions 

_ r+(z,t) _ c 3 Z(t)S(z,t) 
MC) = _ > S (C) = 7 — 77^ . (8a) 



•(*) 



fcr,(t) 



which in the turbulent region z + > z^ uf are assumed 
to be only functions of the "outer variable" £ = = 
z/Z(t). Thus Eq. (4c) can be written as 



K(z,t) = T*(t)k(()/c 2 



(8b) 



We assume also that, the temporal growth of the bound- 
ary layer is proportional to the friction velocity, 



dZ(t)/dt = a\A"*(*), 



(8c) 



with a dimcnsionless constant a. 

We now substitute K(z,t), S(z,t) and r(z,t), ex- 
pressed by Eqs. (8) in terms of k((), s(() and T»(f) into 
the momentum and turbulence energy balance Eqs. (4). 
We also: i) neglect the viscous contribution vS in 
Eq. (4a) in the region z + > z£ n{ ; ii) account only for 
the leading contribution to dK/dt (that originates from 
dZ/dt) and neglect the contribution of dr^/dt. Similarly, 
in Eq. (2) we account only for the leading contribution in 
lnRe r 3> 1 terms. The result is two ordinary differential 
equations for k(Q and s(£) with explicit expression for 
T,(t): 



dk 
d( 
dk 

k— + a(s 
d( 



bk 3 / 2 



j d r~r dk 

d-jjQvk — = bks 
dQ dQ 



.(*) 



ln(i/i* 



(9a) 
(9b) 



In Eq. (9a), the RHS again represent the kinetic energy 
production, and the three terms in the LHS describe tem- 
poral dependence, energy dissipation and the diffusion, 
respectively. The above system of equations for the func- 
tions fc(C) and s(C) indicates that the factorization (8a) is 
consistent with the logarithmic accuracy of the momen- 
tum and turbulence energy balance Eqs. (4) in the sense 
that fc(C) and s(£) depend only on one variable £. 

The boundary conditions for these equations at the 
edge of TBL, Q = 1, are: 



s(l) = 0, fc(l)=0 



(10a) 



To formulate the boundary condition near the wall, it is 
noted that Eqs. (9) are valid only for z + > zt f ~ 50, 
which corresponds to £ > £buf = z^" uf /Re T . In this region 
the mean velocity satisfies the von-Karman law (7) which 
gives 5 + (2^ uf ) = l/(Kz+ uf ). Noting that the full Min- 
Modcl (4) leads in the stationary case to K + (z^ uf ) = c~ 2 
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FIG. 1: Color online. Left: Profile of k((), Eq. (12) with fit constants, presented after Eq. (13) is shown by blue solid line, 
profile fc(C), given by Eq. (9b), - blue dashed line and small £ asymptotics k(C) = 1 — /3oC ~ black dotted line. Profile of £s(£) is 
shown by red dot-dashed line. Middle: Log-plot of mean velocity profile V + (z + ) in the MM for Re T = 610 - blue dashed line, 
Re T = 1150 - red solid line and von-Karman log-law Eq. (7) - thin black solid line. Right: Comparison of the mean velocity 
profiles, given by improved Min-Model (dashed lines) with DNS data (solid lines) for Ree = 1430 (giving Re T sa 610) in blue, 
and Re e = 2900 (Re T w 1150)- in red. 



and to Eq. (7) with k = c 3 jb, we obtain with help of 
Eq. (8a): 



Cbufs(Cbuf) = 1 , fc(Cbuf) — 1 



(10b) 



Since £buf = z£ u{ /Re T , then we take the limit of Eq. (10b) 
as £buf - * for asymptotically large times when In Rc T 3> 
1. It should be noted that boundary conditions (10b) arc 
formulated near the wall (in the log-law region) but not 
at the wall, in the viscous layer, where S + ~ 1. 

With the boundary conditions (10b) and Eqs. (8a) we 
obtain the mean velocity V(z,t) by simple integration: 



V(z,t) = ^(tj 



i r z/z 

f- / s((')dC 



H 



He: 



(11) 



Here the lower limit of integration and term B are chosen 
to satisfy the von-Karman law, Eq. (8c), in the turbulent 
log-law region (for details, see Appendix). 

In order to solve analytically Eq. (8), we now introduce 
a polynomial form of k(C): 

HO = (i - C) 2 (i + 0iC + + AC 3 + AC 4 ) 2 , (12) 

which satisfies the boundary conditions Eq. (10). Sub- 
stituting (12) into (8) and (10), and neglecting terms of 
third-order and higher in the resulting equations give the 
five constants a , 0\ , . . . ft as 

01 = l-0o/2, o = a/K, 

02 = /3o(l + c K /d)/2-(2 + /3 3 +/34), 

3 = 3-/3 (1 + 13ck/V6)/2-2/3 4 , ( i3 ) 

04 = -4 + 0o [1 + c 4 (c 4 + bd)/24 + 371ck/(108cZ)] /2 . 

Substituting k = a J C S (C) <^C [which follows from inte- 
grating Eq. (9b) ] into Eq. (13), the constants 0q , . . . 04 
can be expressed in terms of the model's three parameters 
b, c and d. For the fixed values of b = 0.36, c = 0.53 and 



d = 0.07, the resulting values of 0i are 0o w 0.996, 0\ f» 
0.502, 02 ~ 1.99 , 03 ~ -3,05, 4 « 1.10. 

Figure 1 (left) displays the profile of fc(£) ( blue solid 
line), given by Eq. (12) and the above values of the con- 
stants O . . .04. Substituting this profile in the energy 
balance Eq. (9a) we obtain the profile of (s(Q, shown 
by the red dot-dashed line. Substitution of the resulting 
profile C S (C) m to Eq. (9b) gives the profile of shown 
in Fig. 1 (left) by the blue dashed line. If the initial func- 
tion k(C) gives the exact solution of the balance Eqs. (8) 
then the functions k(Q and k(Q would coincide. Fig. 1 
shows that these functions are reasonably close. Thus, 
we conclude that the simple polynomial form (12) ap- 
proximates the solution of the balance Eqs. (8) in the 
entire TBL, < C < 1, with good accuracy. 

Comparison of Min-Model results with DNS. In 

the present section, we compare the results of our Min- 
Model for TD-TBL with those obtained by DNS of SD- 
TBL [1—3]. This comparison is justified by the fact that 
TD-TBL and SD-TBL are asymptotically equivalent in 
the limit of large Re r [4] . 

a. Mean streamwise velocity. In our solution the 
unsteady and diffusion terms in the turbulence energy 
balance vanish near the wall, for £ — > 0. Accordingly in 
this region the mean streamwise velocity V(z,t) satisfies 
the von-Karman law (7), which gives [in our normaliza- 
tion Eq. (8a)] s(C) 1/C Indeedan Fig. 1 (left), ( s(() 
become constant for small £. 

The figure also shows that the product C S (C) m the 
wide region < £ < 0.6 exceeds unity, which is the level 
of C S (C) m the log-law layer. Accordingly, after integra- 
tion s(() over (, the resulting mean velocity profile V^C), 
shown in Fig. 1 (middle) exceeds the log-law level. Thus, 
the Min-Model clearly demonstrates the wake contribu- 
tion, described by Coles [8] for stationary channel flow. 
The wake contribution is also seen in our DNS results for 
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FIG. 2: Color online. Energy balance in Eq. (9a). Energy production term, k(() s(() - blue solid line, the energy dissipation - 
black dotted lines, the diffusion term red dashed line, temporal term - green dash-dotted lines and the energy input (difference 
between the energy production and dissipation) - black solid line in the insert of left panel. Different panels display different 
region of TBL. 



the spatially- developing TBL. However, the size of the 
wake, given by Min-Model Fig. 1 (Middle) is smaller than 
that in DNS see Fig. 1 (Right). This discrepancy origi- 
nates from the estimate of the outer scale of turbulence £. 
In the Min-Model I is estimated, following von-Karman , 
as the distance to the solid wall z. This assumption is 
valid only for z <C Z, otherwise the scale t is affected by 
the free upper boundary and saturates. Recently [9] the 
effect of saturation of I (at level of 0.3 of the channel half- 
width) on the mean velocity was studied in the channel 
flow. Accounting for the saturation of I in TD-TBL at 
the level ~ Z/2 (which affects only the shear) improves 
the calculated velocity profile, Fig. 1 (Right). The size of 
the wake, given by Improved Min-Model compares well 
with that in DNS. This correction docs not effect the 
other results. 

b. Turbulence kinetic energy. Figure 2 displays in- 
dependence of various terms in the energy balance. One 
can see in the left panel that close to the wall, £ < 0.25, 
the energy production is well balanced by the energy dis- 
sipation; the diffusion and unsteady terms play a minor 
role. They become relevant only for ( > 0.3 as shown in 
the middle panel. Thus the region, £ < 0.25, can be con- 
sidered as "equilibrium" TBL with local spatial energy 
balance. For further clarification, we plot in the insert 
the difference between the energy production and dissi- 
pation. The difference is denoted as the energy input. 
This is the part of the energy flux that is required for the 
temporal development of TBL. Notice that the energy 
input is finite, while the energy production and energy 
dissipation themselves diverge as l/£ for ( — > 0, as shown 
in Fig. 2, Left. 

The energy input is essential only in the first half of 
the TBL, £ < i. As both energy production and en- 
ergy dissipation become smaller (see Fig. 2 Middle and 
Right) their difference vanishes. The energy input can 
be totally neglected in the outer tenth of TBL, where 
C > 0.9. The only source of energy here is the turbu- 
lent diffusion, which transports energy from the region 



C < 0.5 to the region ( > 0.5. Thus, the turbulent dif- 
fusion leads to increase the width of the TBL in time. 
Figure 3, Middle, shows good agreement between the an- 
alytical (dashed lines) and DNS (solid lines) profiles of 
the turbulence kinetic energy. The observed discrepancy 
between the analytical and DNS results for z + < 50 is ex- 
pected since our Min-Model was not designed to predict 
the buffer sub-layer flow. 

c. Reynolds stress is one of the most important 
characteristics of TBL, being responsible for the mechan- 
ical balance in the turbulent log-law and outer regions. 
In stationary regimes, the Reynolds stress is prescribed 
by the outer conditions, maintaining the flow. For 
example, in pressure driven channel flow (of half-width 
L) t(z) oc (1 — z/L), while in the zero pressure-gradient 
plane Couette flow t(z) = const. In developing TBL, 
t(z) is not known a priory, and thus the Reynolds 
stress has to be determined self-consistently. Our model 
[Fig. l(left),fc(C)] predicts that in the first half of TBL, 
near the wall r(z) ps [1 — z/Z], as in the channel 
flow, while in the last quarter of TBL, near the free 
boundary, it decays quadratically, t(z) oc [f — z/Z] 2 . 
This prediction is in a good agreement with the DNS 
data, Fig. 3 (Middle). The production term is also well 
described by our model (see Fig. 3, right panel) 

Summary: We have presented a simple analytical 
model (Min-Model) of the physics of time-developing 
TBL in a Newtonian fluid. The model is based on 
the exact equation for the momentum flux and on 
the model equation for balance of turbulence kinetic 
energy with the production, dissipation and turbulent 
diffusion terms. The Min-Model results in a partial 
differential equation for the kinetic energy, which, in the 
limit of large evolution time was reduced to a relatively 
simple ordinary differential equation. We obtained an 
asymptotic formula for the time dependence of width 
of TBL and approximate analytical solutions for the 
profiles of the mean velocity, the turbulent velocity 
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FIG. 3: Color online. Comparison of analytical predictions (dashed lines) and DNS results (solid lines) for the turbulence 
kinetic energy (Left), Reynolds stress (Middle) and production term (Right) . Blue lines corresponds to Reg=1430 (giving 
Re T S3 610), red lines - Re e =2900 (Re T « 1150). 



fluctuations and the Reynolds stress as a function of 
time and distance from the wall. These profiles are in 
good agreement with the DNS observations. In future 
work we will demonstrate the asymptotic equivalence of 
temporally- and spatially-developing turbulent bound- 
ary layers by the relationship x <^ Voot (here x is the 
distance from the front edge and Voo is the free-stream 
velocity) and will clarify the difference between these 
two regimes in pre-asymptotical region, where lnt + or 
In x + are not very large with respect of unity. 
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Appendix. Here we describe the derivation of Eq. (11). 
From the definition of mean shear (1), the mean velocity 
is written as 



Substituting S(z,t) from (8a) into (14) gives 



V(z, t) = V(z hu{ , t) + f C s(C ') dC ' , (15) 



where £ = z/Z(t) , £buf = zt> u f/Z(t), . The mean veloc- 
ity U(zbuf) at the edge, (zbuf), of the log-law region is 
given by the von-Karman law as 

V(Zbuf) = \ B + In«bufl- ( 16 ) 

Noting that s(Q = 1/C for C < 1, Eq. (15) becomes 



, r I ft 

V(z,t) = B + - s(C')dC 

1 K J Co 



(17) 



V(z, t) = V(z huf , t) + / S(z', t)dz' . (14) where Co = t*/Z(t) = l/Re T 
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